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Abstract Heavy-light QCD currents are matched with HQET currents at two loops and 
leading order in 1/m. A single formula applies to all current matchings. As a by-product, 
a master formula for the two-loop anomalous dimension of the QCD current q^ 1 . . . y^lg 
is obtained, yielding a new result for the tensor current. The dependence of matching 
coefficients on 75 prescriptions is elucidated. Ratios of QCD matrix elements are obtained, 
independently of the three-loop anomalous dimension of HQET currents. The two-loop 
coefficient in /b*//b = 1 ~ 2a s (m,b)/37i — K^al/n 2 + 0(a^, l/m&) is 

K b = — + — 7T 2 + — 7T 2 log 2 - -C(3) - — N t + A C = 6.37 + A c 
12 81 27 S 9 SV ; 54 

with Ni = 4 light flavours, and a correction, A c = 0.18 ± 0.01, that takes account of the 
non-zero ratio m c /mb = 0.28 ± 0.03. Fastest apparent convergence would entail a s (fi) at 
\x = 370 MeV. "Naive non-abelianization" of large-iV; results, via Ni ^ Ni — y, gives 
reasonable approximations to exact two-loop results. All-order results for anomalous 
dimensions and matching coefficients are obtained at large (3q = 11 — |iVj. Consistent 
cancellation between infrared- and ultraviolet-renormalon ambiguities is demonstrated. 
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1 Introduction 



QCD problems with a single heavy quark staying approximately at rest are conveniently 
described by an effective field theory — HQET [1], |2|] (see || for review and references). 
QCD operators are expanded, in powers of 1/m, in terms of HQET operators, m being 
the on-shell heavy-quark mass. In this paper we study, at the two-loop level, the relation 
between currents in the full theory and the effective theory. 

Specifically, we consider heavy-light bilinear currents, J — Qo^Qo, where T is a Dirac 
matrix and the subscript denotes an unrenormalized quantity. The MS-renormalized 
QCD current, J(//) = Zj 1 ^) J , is expanded in HQET operators as 

J( l ,) = Cr(i,)j(ii) + ±J2B i ( f i)d i ( f i)+0(±) , (1.1) 

lit- \ 1 1 L / 

where J(/i) = Zj 1 ( / u) J is the corresponding renormalized HQET current, J = %TQq is 
the unrenormalized HQET current, Qq is a two-component static-quark field, satisfying 
loQo — Qoi an d Oi(fi) are dimension-4 HQET operators, with appropriate quantum 
numbers. The meaning of the operator equality (|1 . 1| ) is that on-shell matrix elements 
of J(/J,), in situations amenable to HQET treatment, after expansion to a given order in 
1/m, coincide with on-shell matrix elements of the right-hand side. 

A single dimension-3 term appears on the right-hand side of ( |1 . 1|) if the Dirac matrix 
r satisfies the conditions 

7 r = ar 7o (<r = ±l), ^V^ = 2ahV, (1.2) 

where h is a function of the space-time dimension, d = 4 — 2e. For an antisymmetrized 
product of n 7-matrices, V = 7^ . . .7^, one obtains h = r](n — 2+e), with 77 = — cr(— l) n . 

It is natural to perform the matching at a scale fi ~ m, where the matching coefficient 
Cr(/i) contains no large logarithm. One can then use the renormalization group to relate 
QCD and HQET currents renormalized at arbitrary scales, fi and fl: 



cxp 



\ a s (m) ) \ a s {m) 



where 7j = d log Zj/d log \i and 7j = dlogZj/dlog/x are the QCD and HQET current 
anomalous dimensions, and (3 = — dloga s /dlog/i. 

To calculate the matching coefficient C*r(/i), we consider an on-shell QCD matrix 

/ \ 1 / 2 1 

element, M{ji) = {ZqZ° s \ Zj (fi)T , obtained from the bare proper-vertex function, 
To, by renormalizing the current and performing on-shell wave-function renormalization. 

To zeroth order in 1/m, it is equal to Cr(/i)M(/i), where M(fi) = ^ZqZ° s ^J ^ Zj 1 ( / u)f . 
Hence we obtain 

Cr( ,).(s) 1/2 |44&, ( m) 

\z%) Zj(p)r ' 

with a fi dependence coming only from Zj{\x) and Zj([a), in agreement with ( |1.3| ). 
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The matching coefficient Cr(/x), for any T of the form ( |1.2[ ), was obtained at the 
one-loop level by Eichten and Hill IJ. The 1/m suppressed matching coefficients Bi(fi) 
in (|1 . 1| ) were obtained for vector and axial currents, at one loop, in [|], ||. Here we shall 
find the two-loop correction to the leading matching coefficient Cr(/i). It is required for 
a more accurate extraction of QCD matrix elements, such as /b and /b*> from HQET 
results obtained, for example, from lattice simulations or sum rules. 



As can be seen from (|1.3| ), it is in general necessary to use three-loop anomalous di- 
mensions, 7j and 7j, in conjunction with a two-loop matching coefficient, Cr{m). Whilst 
some QCD anomalous dimensions are known at the three-loop level, the universal HQET 
current anomalous dimension, jj, is currently known only at the one- || [7j and two- 
loop H, 01 [T(J| levels. There are, however, three cogent reasons for calculating Cr(m) to 



two loops. First, previous experience O, 12] strongly suggests that two-loop finite effects 



dominate, numerically, over three-loop anomalous dimensions, near a heavy-quark mass- 
shell. Secondly, the HQET anomalous dimension is independent of the Dirac structure T 
and hence is absent from important ratios of physical quantities, such as /b* //b- Finally, 
the three-loop term in jj may be calculated in future. 

Motivated by these considerations, we here compute Cr(m), to two loops, for an 



arbitrary Dirac structure of the form (|1.2|) . Section 2 gives our method and general result, 
for any T. In Section 3 we consider ratios of matching coefficients. After elucidating the 
dependence on 75 prescriptions, we show that the two-loop correction to /b* //b, at zeroth 
order in 1/m, is comparable in size to each of the effects previously computed, namely 
the one-loop correction of [p]] and the 0(l/m) correction of [[13], |H| . In Section 4, we 
calculate anomalous dimensions and matching coefficients to all orders in a s , in the limit 
of a large number of flavours, Nf. We extend recent analyses [|l|, [16, [L7| by obtaining 



the MS matching coefficient for an arbitrary T in this limit, confirming the cancellation 



of renormalon ambiguities in physical matrix elements |17| and validating an additional 



consistency condition. Section 5 gives a summary of our main results. 



2 Two— loop matching calculation 

The bare proper vertices, Tq and To in (|1.4j ), may be evaluated at any on-shell momenta 
of the light and heavy quarks. The calculation is greatly simplified by choosing the mo- 
mentum of each quark to vanish in HQET. This corresponds to a heavy-quark momentum 
p = mv in QCD, where v = (1,0) in the heavy-quark rest frame. One- and two-loop 
diagrams for the vertex are presented in Fig. [I]. 

In the quark loop of diagram [l|b we include Nf = Ni + 1 quark flavours, one of which 
is the external heavy flavour, whilst the remaining Ni light flavours of quark have masses 
rrii < m, any of which may be zero. In the case of an external b quark, Nf = 4 + 1 = 5 
and only the mass ratio m c /mfc ~ 0.3 is significantly greater than zero. 

We stress that the loop containing the external heavy-quark flavour, with mass = 
m, must be included in all 6 terms on the right-hand side of ( |1.4| ). This loop is certainly 
present in the QCD terms and we must therefore include it the HQET terms, as well. 
The HQET modification of QCD propagators occurs only along a single heavy-quark line, 
going through all diagrams; HQET is blind to the contents of loops. In effect, HQET has 
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no memory that the static-quark line had one of the massive flavours that still occur 
in loop corrections. Just as in conventional QCD, heavy loops decouple from physical 
quantities, at low momenta; not from renormalization constants, or bare vertices. We 
shall demonstrate the physical decoupling of the t-quark loop from / B * / /b and show that 
the contribution of the 5-quark loop is considerably less than that of a light-quark loop. 

We proceed to analyze, in turn, the 3 HQET quantities and 3 QCD quantities on the 
right-hand side of ( |1.4|) . 



2.1 HQET calculation 

Here we deal with the HQET quantities Tq, Zj, and Zq. 

In the case of the HQET proper vertex, T , all the dimensionally regularized diagrams 
of Fig. |T| vanish, for external quarks with zero residual momenta. This is immediately 
obvious for all except diagram [l|b, with a massive-quark loop, since no other diagram 
contains a scale. In fact, this diagram vanishes too, because the gluon propagator contains 
a quark-loop insertion and is hence transverse, i. e. proportional to (g^ — k^k v /k 2 ). This 
tensor is contracted with v u , from the heavy-quark vertex, and with ft'y fJ- /k 2 , from the 
light-quark vertex, 7 M , adjacent to an internal light-quark propagator, ft/k 2 , with the 
same momentum as the gluon. The resulting contraction gives — 7 ■ k'-fo/k 2 , which is odd 
under reflection of the spatial momentum k and hence vanishes when integrated over the 
loop momentum k. Thus T = Tb, where Tb = u q TuQ is merely the Born term. 

The HQET current renormalization Zj does not depend on V. In the MS scheme, at 
two loops, it is fully determined by the anomalous dimension || [| |i~0f : 



7j 



2 r 



Cr ( 2 " 16 C(2) 



49 
~6~ 



4C(2) 



+ ^T F N f 



(2.1) 



where Ca = 3, Cp = 4/3, Tp = 1/2, in the case of QCD, and we omit terms, here and 
subsequently, without further comment. 

The HQET on-shell wave-function renormalization constant is defined by Zq = 

(l - (dEo/do;^ =o) where S (^) is the bare static-quark self-energy, at residual en- 
ergy ijj. As in the case of the vertex, only diagrams containing a loop with a quark of 
non-zero mass can contribute, since no other diagram contains a scale. Thus we obtain 



dS 



duj 



-iCpgl 



oj=0 



d d k v»v u n^O) 
{27c) d {vk) 2 k A 



—i 



C F g 2 f d d k 



k 2 



d-lJ {2n) d \{vky 



ng(fc) 
k e 



(2.2) 

where li^k) = (g^ u — k^k u / k 2 )IY^{k) /(d—1) gives the contribution of quark loops to the 
polarization operator, whose trace depends only on k 2 . We now average over all directions 
of k, in (i-dimensional space, obtaining < k 2 / (vk) 2 > = 2 — d. A dimensionally regularized 
result is then readily obtained from two-loop massive bubble integrals [DJ, which give 



dSr 



du 



— 8C pTp 



uj=0 



\^ n J m,i>0 



(d-l)(d-6)T 2 (l + e) 
(d - 2){d - 4) 2 {d - 5)(d - 7) 



(2-3) 
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whose Laurent expansion yields the on-shell wave-function renormalization constant 



m 



16 , rrii . 26 
^ + -log^ + C(2) + - 
ji 6 n 9 



(2.4) 

We note that (|2.4|) exhibits non-uniformity as wij — > 0: quarks with masses m 8 > are to 
be included; any quark with rrii = is excluded, by dimensional regularization. This non- 
uniformity has the same infrared origin as that observed in on-shell QCD wave-function 
renormalization ||12||, to which we now turn. 



2.2 QCD calculation 

For continuity, we deal with the QCD quantities, To, Zj, and Zq, in reverse order. 

A gauge-invariant on-shell wave-function renormalization constant, Zq, was obtained 



in d dimensions, at the two-loop level, in [0. A quark flavour with a small but non- 
zero mass gives a contribution to Zq different from that of a zero-mass flavour. The 
same is seen to be true for Zq in ( |2.4j ). Each discontinuity originates from the infrared 
region of small gluon momenta, where HQET does not differ from QCD. Therefore these 
discontinuities should cancel in the ratio Zq / ' Zq, which we indeed find to have a smooth 
limit, as rrii —* 0. This uniformity provides a strong check of the massless and massive 



quark-loop contributions to Zq, obtained in [fL2j| . The ratio 



ZqIZq l + 2C F T F ^) Y.^ -) (2-5) 



(^Q/^Q)rrii=0 V47T/ i=1 



III 



may be expressed as an integral of the difference, II (—mf/k 2 ), between the polarization 
operator with a quark of mass rrii in the loop and that with a massless quark. This 
difference contains neither ultraviolet nor infrared divergences. The integral over the 
gluon momentum k is likewise well-behaved, and may be performed in 4 dimensions. As 
is the case for all two-scale calculations in this paper, the result may be expressed as a 
combination of three basic dilogarithmic integrals, defined and evaluated in Appendix A. 
Combining Zq, from | 12fl , with Zq, from ( |2.4|) , we obtain 



A z (r) = -4Ax(r) - 2A 2 (r) - 12A 3 (r) . (2.6) 

The QCD current renormalization constants Zj are known to two loops for some, but 
not all, of the currents that we study. We shall simply apply a minimal Ansatz for Zj and 
extract its coefficients from the requirement of finiteness of (O), checking our resultant 
master formula for jj against known special cases. 

The most difficult part of the problem is the calculation of the bare QCD proper 
vertex, r . We shall first find it in the case when all light flavours are massless, and 
include later the effects of finite light-quark masses. Since we calculate the bare vertex 
on the renormalized mass-shell, p = mv, it is convenient to express the result in terms of 
the on-shell mass m, treating the one-loop G?-dimensional counterterm, Am = m — mo, 
as a perturbation. 
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Initially, we make no assumption about the properties of the matrix T, and note that 
each two-loop diagram for Tq may be written as a sum of terms, each of the form 

where / is some integral over loop momenta, / is even, and l+r < 8. After the integration, 
jfi 1 ...fii-u 1 ...u r can con t a j n on ly giw anc [ v a _ Resulting contractions of pairs of 7-matrices 

on the left, and of pairs on the right, merely produce additional terms of the form ( |2.7|) , 
with smaller values of I + r. Before performing the remaining contractions, one may 
anticommute 7-matrices, so as to arrange that ^ occurs only on the extreme left, or the 
extreme right, with the contracted indices in between occurring in opposite orders on 
the left and right of 17 Additional terms, arising from anticommutators, have fewer 7- 
matrices, with / remaining even. Repeating this procedure for all values of l+r, from 8 
down to 0, we may cast any diagram in the form 

D = u q [V(x + xif) + ^7^7° (x 2 + x 3 tf) + 7« 7/^7^7° O4 + x$i>) 

+ ^7a7/37 7 r7 7 7 /3 7 a (^6 + x-jj) + 7a7/3777<5r7 <5 7 7 7 /3 7 Q 2 ; 8]«Q • (2-8) 



Now we assume that the matrix Y has the properties (|1.2j ). The effect of each con- 
traction is then to produce a factor 2oh. Terms with an odd number of contractions 
necessarily contain ji on the left, which yields an extra a when moved to the right, where 
it merely gives $uq = uq. Thus we obtain a result involving only powers of h: 

D = [{x + zi) + [x 2 + x 3 ){2h) + (x 4 + x 5 )(2h) 2 + (x 6 + x 7 )(2hf + x 8 (2h) A ] T B . (2.9) 

We can find the coefficients x + x\, x 2 + x 3 , x 4 + x 5 , x 6 + x 7 , x$, for each diagram, by 
taking, separately, a trace of the 7-matrices on the light-quark line with Li, and a trace 
on the heavy-quark line with Hi, using the following 5 forms of Li x Hf 1 x (1 + jf); 
7^ x (1+^)7^; 7^ x (1 + ^)7^; liflv^pi x (l + ^)7 p 7 l/ 7 /i ; j^vlpla x (1 +^h a YYl fl - 
Performing the same operation on the generic form ( |2.8| ), we obtain 5 equations relating 
these double-traces to the desired coefficients. Inverting these equations, once and for 
all, we may then convert the integrand of any diagram into a 4th order polynomial in h, 
with coefficients that are scalar functions of v and the two loop momenta. For reliability, 
we checked our general result, for each two-loop diagram, against brute-force evaluation 
of 8 specific cases of Dirac matrix. The general one-loop diagram |l]a, and its associated 
mass-counterterm contribution, were similarly reduced to quadratic functions of h. 

All the resulting one-loop scalar integrals can be expressed in terms of 

i r d d k r(l + e) 



'0 



7T" 



d l 2 J k 2 + 2vk e(l - e) ' 



(2.10) 



The two-loop integrals may be reduced to combinations of Iq and two further terms: 

1 r d d kd d l 
h — — 1 



7C d J k 2 {l- k) 2 {l 2 + 2vl) 

1 - 4e r(l + e)r 2 (l - e)r(l + 2e)r(l - 4e) 

~ 2£ 2 (l-2e)(l-§£)(l-3e) r(l-2e)r(l-3e) ' 

/=_!/■ ^ (oil) 

2 ir d J (k 2 + 2vk){l 2 + 2vl)((k + l) 2 + 2v(k + l)) K ' ; 

3(5d- 18) (d- 2) 2 2 2(d-4) 16(d-4) 2 

-Jo 771 z-Ji - -nr. 3T7o3 T7a 1 \ £ ) > 



2(3d- 8)(3d- 10)(d- 3) u 2d - 7 (3d - 8) (3d - 10)' 
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where [1181, IIU 



1(e) = / + 0(e), / = 7r 2 log2-|C(3). (2.12) 



To achieve this reduction, we used the package RECURSOR p9| , written in REDUCE 
to implement recurrence relations derived from integration by parts |fLl| , [12| . 

So as to have a strong check, we performed the calculation in an arbitrary covariant 
gauge, verifying that the sum of all diagrams for the on-shell vertex, including the mass 
counterterm, is gauge invariant in d dimensions. The result can be written as 



" " Cf WW W^) + F W h hh jik 1 j k ' ( } 



in terms of 4 colour factors, 3 integral structures, and 4 current-specific functions. The 
colour factors are chosen as follows: Co = Cf — \Ca\ C\ = Cf] C*2 = TpNi, from iVj 
massless loops in diagram |l]b; C3 = Tp, from the heavy-quark loop. The integrals 

j = ( d - (2 u) 

8{d- l)(d-3) 2 (d-4) 2 (rf-5)(d-6) ' l ' ; 

J l — TTTi TT71 7W71 7T77T1 5 J 2 



}(d-l)(d-3)(d-4)(2d-7) ' 4(d - l)(d - 4) 2 (rf - 6) ' 

are chosen to make the coefficients ajik polynomials in d. All dependence on T resides in 
hk = h k , for k = 0,1, 2, and = h 3 (h + d — 4). Only diagrams |l|f and |l]i have sufficient 
7-matrices on each side of the vertex to generate h 3 and /i 4 terms, and we find that each 
diagram produces the combination h 3 , for each integral structure. Moreover, the integral 
J2 arises only from diagram |T]j, with colour factor Co, or from the heavy-quark loop in 
diagram |l|b, with colour factor C3. Thus several of the 48 coefficients ajik vanish. The 28 
that survive are listed in Appendix B. 

We complete the calculation by including the effect of non-zero light-quark masses. 
The difference of diagram |l]b, with a quark of mass m; in the loop, and the same diagram 
with a massless quark, contains neither ultraviolet nor infrared divergences. Therefore we 



calculate it in 4 dimensions, obtaining a further term, to be added to (|2.13|) , of the form 

T7 = c ' r '(s)£>Ur). (2 ' 15) 



where the sum runs over light flavours only, and Ap can be reduced to the integrals ([A. 2 ) 
as follows: 



Arfrl 



4 r 



2(l-h 2 )A 1 (r)-(l + 2h-2h 2 )A 2 (r) + 2h(2-h)A 3 (r) , (2.16) 



with Ap(0) = 0, by virtue of its definition. As a strong check on the coefficients of the 
colour factors C3 and C2 in ( j2.13j ), coming from loops with quarks of masses rrii = m and 
rrii = 0, respectively, we have verified that setting d — 4 in their difference agrees with 
setting r = 1 in ( p,16|) . 
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2.3 Results 



We re-express the bare QCD vertex (p. 13 ) in terms of the MS coupling, ct s (/i), using 
one-loop coupling-constant renormalization, with Nf = Ni + 1 flavours, and expand the 

and 



result as a Laurent series in e. Using Zq from (T^|, T = r B , Zj(fj) from (|2J 
Zq from fl2.4| ), we obtain the generic QCD current renormalization constant, Zj(fi), by 
requiring the finiteness of (|1.4|) . The corresponding anomalous dimension is that of the 
QCD current J(fj) = qj^ 1 . . .j^q, which does not depend upon whether the quarks are 
heavy or light. We find that 



7J = _<7 *(n-l)(n-3) 
Zty 



1 + ^(c/^L 



19 



C A 



3(n-2) 



19 



-^(£) 2 <" 



l)(n- 15) 



11C A - ATpN 



f 
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(2.17) 



with a check being provided by the absence of the signature rj, which merely distinguishes 



different components of J. We reproduce the known result for the scalar current pi 
with n — 0, and, of course, obtain a vanishing result for the vector current, with n — 1. 
The two-loop anomalous dimension of the n = 2 current, qcr^ul, appears to be absent 
from the literature, though it could be derived by straightforward calculation in massless 
QCD. Instead, we have obtained 



U\n=2 



4?r V 2 18 9 } 



(2.18) 



as a by-product of a more difficult massive on-shell calculation. The currents with n = 3 
and n = 4 are non-singlet axial and pseudoscalar currents, with the 't Hooft-Veltman 75, 
whose anomalous dimensions are known to differ, at two loops |22, 23 and beyond |2^, |24]|, 
from those with the naively anticommuting 75, with n = 1 and n = 0, respectively. This 
difference is apparent in the final term of ( [2.17| ), which lacks invariance under n — > 4 — n. 
Our results for n = 3 and n = 4 agree with E2] and EB[, respectively. 



Finally we arrive at the general expression for the two-loop matching coefficient 



C v {m) = l + C F 



47T 



3(n-2) 2 + (2-rj)(n-2) -4 



Cp(ip + C A a A + T F y2<a f + Aj ( — 



(2.19) 



where the sum now includes the heavy flavour, with = m. The one-loop term coincides 
with The coefficients in the two-loop term are 



a F = - yC(2)) (n - 2) 4 + ll(n - 2) 3 - ^ v (n - 2f 



+ 



~ + 48C(2) - ^/ ) (n - 2) 2 - 2?7(n - 2) 2 - 20(n - 2) 
o o 



+ ( y - y C(2) + \l) v(n-2) + 6 -^- 81((2) - 8C(3) + 12/ 
43 4 \ 

'M [ —— + -C(2)J (n - 2) 4 - 2(n - 2) 3 + r/(n - 2) 3 



(2.20) 
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+ 



V216-f^ + 3 J j^- 2 ) 2 
~ + 8C(2) - V (n - 2) - 



143 

" ~18 
29017 

432 



(n-2) 



+ 29C(2) + 2C(3)-6J 



a f 



445 8 
3 



54 



C(2)) (n - 2) 



> 2, nN 38 . 1745 20 _ 

2 -(n - 2) + — r?(n - 2) + 1 C (2) 

g v ; -r g /v JT 108 3 SU 



The mass correction Aj(r) = Ap(r) + A^(r) depends only on /i| e= o = 7/(ti — 2), be- 
cause ( |2.6|) and Q2.16j ) result from finite, 4-dimensional integrals. Its general form is 



Ar(r) 



3 

2 (2(n 



2 (2(n - 2) 2 - 5) Ai(r) + (4(71 - 2) 2 - 477(71 - 2) - 5) A 2 (r) 



4r/(n - 2) + 9) A 3 (r 



:2.21) 



The heavy-flavour contribution, with m 8 = m, is proportional to a/ + Aj(l), where 



Aj(l) = - [2(7 + 2C(2))(n - 2) 2 + 4(-5 + 2C(2))r/(n - 2) + 17 - 38C(2) 



(2.22) 



is always opposite in sign to af and close to it in magnitude, resulting in a heavy-flavour 
contribution that is always small, in comparison with that from a light flavour. 

Table 1: Matching coefficients 



r 






Cr 


m) 






1 


1 + 


2 a s (m) 

3 7T 


+ (10.92- 


O.6ON1 + 0.04) ( 


v T , 




7o 


1- 


2 cc s (m) 

3 7T 


-(4.20- 


0.44iVj + 0.07) ( 


/ a g . N 


f 


7i 


1- 


4 «s(m) 

3 7T 


-(11.50- 


0.79iV, + 0.09) ( 


o s 

V T , 




7o7i> 7i72 


1- 


4 a s (m) 

3 7T 


-(16.19- 


1.13iV| + 0.13){ 


Os 




7o7i72 


1 




-(10.98- 


0.77^ + 0.11) ( 


/ a g . N 


f 


7i7273 


1+ 


2 a s (m) 

3 7T 


-(2.78- 


0A2N t + 0.08) ( 


'as' 
v T , 




7o7i7273 


1 + 


10 a s (m) 

3 7T 


+ (19.75 


0.64JVi + 0.00)( 


/ a 5 . N 


f 



There are 8 distinct matrices, T, in 4-dimensional space-time: antisymmetrized prod- 
ucts of up to 3 spatial 7-matrices, and the same products multiplied by 70. For each 
current, we obtain a matching coefficient Cr(m), by choosing appropriate values of n and 
77 in the formulae (|2.19 ) to ( 2.22j) . With N\ zero-mass quarks, and QCD colour factors, we 
obtain the numerical values of Table [T], two of which coincide (for reasons explained in the 
next section). Typical matrices V are presented in the first column; Cr(m) is, of course, 
the same for other spatial components. Coefficients of (a s /7r) 2 in the two-loop corrections 
are written as sums of three terms: the quenched contribution, without quark loops; the 
light-quark loop contribution; and the heavy-quark loop contribution. One can see that 
the quenched contribution is usually of order 10; each light flavour gives a contribution of 
order 1, with the opposite sign; and the external flavour contributes of order 0.1. In most 
cases, the two-loop correction has the same sign as the one-loop correction. 
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3 Ratios of matching coefficients 



First we consider ratios that reflect differences in 75 prescriptions. Then we analyze ratios 
of meson matrix elements. 



3.1 't Hooft— Veltman versus anticommuting 75 

It is generally accepted that there exists no (i-dimensional generalization of 75 that an- 
ticommutes with all 7^. Nor is there any covariant (i-dimensional generalization of the 
Levi-Civita tensor, e^ uap . Instead, one may use the 't Hooft-Veltman 75, 

7 HV = i7 7 l 7V; {31) 

which makes multi-loop calculation rather time-consuming, even when implemented co- 
variantly, using antisymmetrized products p3| p4| , p5| p6| . However, there is also a general 



belief that one may use a naively anticommuting matrix, 7^ c , in open fermion lines, and 
in loops with an even number of 75-matrices, without encountering contradictions. 

The antisymmetrized products of n and 4 — n 7-matrices are related to each other 
by multiplication by 7^ v , and have different matching coefficients. On the other hand, 
multiplying V by the naively anticommuting 7^° does not change h in (|1.2| ), and hence 
cannot change the matching coefficients. Therefore ratios of matching coefficients with n 
and 4 — n antisymmetrized 7-matrices express differences between using the different 75 
prescriptions. 

HQET matrix elements do not depend on the 75 prescription. For example, 

<0|g 75 HV Q|B> = <0|g 75 AC Q|B> . (3.2) 



To prove this equality, consider the HQET operator-product expansion p7j f° r the cor- 
relator of two pseudoscalar currents. In every contributing diagram, the static-quark 
propagator between the two 75-vertices has the 7-matrix structure (1 + 7o)/2, which is 
unaffected by gluon vertices. Since every sort of 75-matrix anticommutes with 70, we can 
move one 75-matrix from its vertex and annihilate it with the other, leaving (1 — 7o)/2 
on the heavy-quark line, independently of the prescription. Since the correlator is inde- 
pendent of the prescription, so is its spectral density, and in particular the ground-state 
contribution. Hence we arrive at (13T 



In contrast to this, the QCD matrix elements do not coincide. We have already 
remarked that the MS-renormalized pseudoscalar currents, J HV (/x) = Z^ilA^lf^ Qo 
and J AC {jj) = Z^Q(/i)go7^ c Qo, have anomalous dimensions that differ, at two loops and 
beyond. The currents are related to each other by a finite renormalization, J AC (/i) = 
Zp(fi) J HV (/x), whose term of order may be found either by comparing renormalized 
matrix elements, at L loops, or more demandingly, by equating d log Zp/d log to the 
difference, 7jhv — 7jac, of anomalous dimensions, at L + 1 loops. The consistency of 
these methods has been demonstrated, at L — 2, by evaluation of massless three-loop 



diagrams | 23| , which yield 



ZM = 1 _ 2Cf ^i + Cf ( S ) 2 m ±pm + 0{a!) , p.. 
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using either method, with a further 0(«g) term determined solely by the first method. 
Consistency demands that this result be obtainable by comparing matrix elements with 
quarks of any mass and momenta, and in particular that 

<0|J AC Qu)|B> CM 
P[N <0|JHV (;U )|B> C 707l7273 (/i) ' [6 - } 

given the equality (|3.2| ). Setting // = m and using (|2.19|) , we indeed verify ( |3.3|) . In 
particular, the two-loop finite-mass corrections (|2.21| ) cancel in the ratio ( |3.4|) , as should 
occur in a result obtainable [23] from three-loop mass-independent anomalous dimensions. 

Similarly, the MS-renormalized 't Hooft-Veltman and naively anticommuting non- 
singlet axial currents are related by a finite renormalization, J AC (/u) = J^ v (/i), 
which has been computed at the one- f2~2|| , two- p5| , and three-loop [^3|, Q levels. We 



are informed by S. A. Larin that the discrepancy between |23| , |2J] and p5| at two loops 



is attributable to the use of the so-called G-scheme in lESfl. The MS result is [E3 



7T \An J 9 

which we verified by calculation of two ratios 

M = C (a) = C (a) ' (3 " 6) 

^ 71 7273 vAv ^ 70 71 72 Ifv 

Note that the weak axial current is J AC . Only for this current does QCD renormalization 
preserve a V — A structure. Thus measurable matrix elements, such as /b, are obtained 
from J AC , not from the chiral-symmetry-breaking 't Hooft-Veltman counterpart. 

The QCD tensor current, = qc^q, is a special case, because inclusion of 7^ v , to 
form J^f = q r ^<J^ v q 1 is merely a space-time transformation, in 4 dimensions, giving, for 
example, J^ v = — U23. Thus there can be no difference of anomalous dimensions between 
J^ v and Jf^v ■> an d hence no non-trivial finite renormalization, relating the two 

prescriptions. Equivalently, we have two methods of obtaining Z T (/u,): 

C7273 (A*) ^7071 (aO 



which proves both that Z T (fi) = 1 and that C 7o7l (m) = C 7273 (m), as found from ( [2.19| ) 



It is easy to see how this equality comes about, purely within our calculation of matching 
coefficients. With n = 2, we have h = rje, so that the only possible origin of a difference 
between C 7o7l (/i) and C 7273 (/i) would be a term h/e in r . Such a term would give an 
77-dependent QCD anomalous dimension, for n 7^ 2, and hence cannot occur. 



3.2 Ratios of meson matrix elements 

After eliminating the currents containing the 't Hooft-Veltman 75, we are left with 4 es- 
sentially different currents. In this subsection, we discuss the ground-state-meson matrix 
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elements: 



<0|(<?75Q) M |B> = -im B f£(ji) , 
<0|g7«75Q|B> = i f B p a , 
<0|g7aQ|B*> = im B *f B *e a , 
<0|(g(x a/3 (5) A1 |B*> = f B *{ij){p a ep -ppe a ) 



(3.1 



where a^ v = |[7 M , y„], and we use, from now on, the naively anticommuting 75 that is 
appropriate to matrix elements of non-singlet weak currents. The corresponding HQET 
quantities coincide: /g = f B = f B * = /b*> as consequences of the heavy-quark symmetry. 
Hence, ratios of QCD matrix elements are equal to ratios of matching coefficients. These 
ratios do not depend on the unknown three-loop HQET anomalous dimension. Similar 
formulae hold for P-wave + , 1 + mesons, if one inserts an extra 75 into all currents. 

It follows from the equations of motion that 

/b(aO = <0|(g75QV|B> = _mB_ ^ g ^ 

/b <0|g7 5 7oQ|B> m(/i) ' 

where is the MS running heavy-quark mass. To leading order in 1/m, we may 

replace m B by the on-shell mass m, obtaining 



f B (m) m Ci(m) 



/b 



m{m) 
2 



C 70 (m) 



1 + 



a s (m) 



7T 



Ct. 



121 



+ 30C(2) - 8/ + C A 



1111 



24 



8C(2) + 4J 



r?) 



(3.10) 



1 + 



AacJm) 



+ (16.01 



7T 



1.04A/i + 0.10) f 



where the final numerical result shows, as in Table |l|, the contributions of: the quenched 
term; iVj massless-quark loops; and the heavy-quark loop. We confirm the analytical 
result of fLl| , in the case of iVj = iV/ — 1 massless quarks, and correct the omission of a 
factor Cf = 4/3 from the numerically much smaller effects of finite light-quark masses in 
Equation (17) of |TT . 



Our most important result is for the ratio of two observable matrix elements: 

<0|g7Q|B*> m B *f B *e 



(3.11) 



<0|<?7o75Q|B> m B f B 
which is independent of the renormalization scale. At leading order in 1/m, we obtain 



/b* 
/b 



C 70 (m) 



An) 



= 1-0! 

1 r 1 
C F -(31 



2tt 

128C(2) 



16/) + C A -(- 



•263 + 144C(2) - 247) 



16^ ^ 

° i=i 



3 7T 



19 
12 



A 9 



m 8 
m 



2 A, 



m 
m 



(3.12) 



(7.75 



0.35iV, + 0.03) ( 
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Note that this is numerically very close to (m(m)/m) 1 ^ 2 , as can be seen from ( |3.10| ). 



At any finite order of perturbation theory, the ratio ( |3.12| ) is determined by a combina- 
tion of loop integrals that contains neither infrared nor ultraviolet divergences. Therefore 
all loop momenta are of order m. In order to demonstrate the absence of infrared sen- 
sitivity (which may sometimes be overlooked in dimensional regularization) we repeated 
the one-loop calculation with a gluon mass, A, finding that the only effect is to multiply 
the 0(a s ) term in ( [Hp by 




OO i-fO / i i \ ^-v-> ^ * 2 



1 r dk 2 „ f k 2 \ . 2vr A ^ / X 2 , 

— - + — , (3.13) 



3 J k 2 + A 2 \ m 2 I 3 m \ m 2 

F(x) = (x + l)yJx(x + 4) - x(x + 3) 



2y/x, x C 1 
2/x, x > 1 



In general, we find that all ratios of matching coefficients are infrared-safe, though those 
involving QCD currents with anomalous dimensions clearly require ultraviolet regulariza- 
tion. 

We now demonstrate that the existence of a very heavy flavour, such as top, has no 
effect on the observable ratio /b* //b- To prove this decoupling theorem, we rewrite ( |3.12|) 

as 



/b* , „ af f] (fx) , _ fa 



2 



W (//) = — CU log - - — 2> 2 - - 7^ + A 2 H - 2 A 3 — 



3 jj 3 ~[ { /i 12 \m/ \m 

where A is a constant, specifying the quenched term in ( |3.12| ), and we have transformed 
the MS coupling to an arbitrary scale /z, using the one-loop iVy-flavour /3-function. The 
explicit //-dependence in fl3.14| ) at 0(a 2 ) is cancelled by the implicit dependence at 0(at s ). 



Now, suppose that we take account of the existence of a super-heavy quark, with mass 
rrih ^> m. The effect is merely to replace Nf by Nf + 1 in ( |3.14p , where 



2?r 2tt 
ensures that a^ f (m/J = «s (m^), and 




(3.15) 



flW+lfc.) = B™ M - 15r F {log™ - H + A 2 - 2A 3 (3.16) 

3 [ /i 12 \m/ Vm/J 



includes the loop-effect of the super-heavy quark in diagram [l|b. Referring to (|A.5| ), we 
find that A 2 (r) — 2A 3 (r) = logr + 1| + 0(r _2 logr), for r = rrtfjm ^> 1. Hence the 



terms in braces in ([3.15D and (|3.16 ) cancel, at any scale \x. This also demonstrates that 



one must match the MS couplings at /i = rrih, as in ( |3.15|) , and not at, say, /i = 2m^. 
Interestingly, the external heavy flavour is also numerically unimportant. Its contribution 
to ( |3.12| ), relative to that of a massless quark, is suppressed by a factor 

1 - g {A 2 (l) - 2A 3 (1)} = U^ 2 ~ 41) = -0.08 . (3.17) 
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Comparable suppressions occur in all our results. 

Turning to the numerical significance of our two-loop result, we substitute Ni 
in (|3.12|) , and add the 1/m correction of | 



14 , obtaining 



/b* 
h 



2 at s (m b ) 
3 



7T 



(Is 



7T 



2 1 

(6.37 + A c ) + — 



A - 8G 2 {m b ) ) + 0(«s, a s /m fe , 1/m 2 ) , 

(3.18) 
0.18 ±0.01, 

for a mass ratio m c jm b = 0.28 ± 0.03, obtained from rrib = (4.8 ± 0.2) GeV, using 
mb,c + A + (/i 2 — imq) / (2rrib,c) and me* = m-Q ± 2fi 2 G / {?>mb) with a kinetic term 



with a finite c-quark mass correction A c = |(A2 {m c /mb) — 2A3(m c /mfc)) 



Hi = (0.5 ± 0.1) GeV 2 . The current world average a s {m z ) = 0.117 ± 0.005 g§, evolved 
down to m& using the three-loop formula, gives a s (rrib) = 0.215 ± 0.018. Therefore 
the one-loop correction is —(4.6 ± 0.4)%, and the two-loop term is comparable to it: 
— (3.1 ±0.5)%. Another way to state the slow convergence of the perturbation series is to 
say that the fastest-apparent-convergence scheme, with a vanishing two-loop correction 
in Q3.14D , would require one to evaluate a s (/i) at far too low a scale: fi = nib/ 13 = 
370 MeV. 

The chromomagnetic interaction matrix element G^^ft) is not well known; sum-rule 
estimates range between G 2 (m 6 ) = —26 MeV \ T3j , and G 2 (^&) = ±21 MeV [Q, with 
large uncertainties. Including the meson residual energy, A, one obtains a 1/m correction 
of ±(11 ± 3)%, from |3], and ±(4 ± 3 ± 2 ± 2)%, from 0, whose three sources of 
uncertainty arise from: sum-rule fitting; an unknown two-loop anomalous dimension; 
and an uncertainty in A. Combining our two-loop radiative corrections with this range of 
1/m effects, we arrive at /b*//b = 1.00 ±0.04, with an uncertainty reflecting the difference 



between |13[ and |14[]. It is, however, quite unclear whether the perturbation series in a s 
is converging fast enough for us to neglect the unknown three-loop term. 



For the last ratio, 



<Q\(qiaoQ)^> 

<0|g 7l Q|B> 



(3.19) 



we obtain 



/b. 



C 707l (m) 
CL (m) 



c F 



= i 

1 /307 
3 VlT 



70C(2) ± 8/ 



C A - 



4277 



± 24C(2) - 4/) 



205 
144 



C(2)±A: 



mi 
m 

2 



A, 



m, 



m 



(3.20) 



1 - (4.69 - 0.34V, ± 0.04) 



The one-loop term vanishes at /x = m, but will of course appear at any other fx. If one 
wants /b*(aO at a scale fi widely separated from m, one needs to use the QCD three-loop 
anomalous dimension of the current qa^q, in conjunction with ( |3.20| ). This anomalous 
dimension is not known at present, though it could be calculated using standard massless- 
quark methods. 
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Finally we note an intriguing pattern in the results ( |3.10| ), ( |3.12| ) and ( |3.20| ): in each 
a reasonable estimate of the two-loop term is obtained by the simple device of replacing 
Ni by Ni — y in the easily computed light-quark-loop contributions. We call this device 
"naive non-abelianization" , since it is based on the hope that results in the non-abelian 
theory may be estimated merely by replacing the leading term in the abelian large- Nf (3- 
function by its non-abelian counterpart. Comparing naively non-abelianized estimates 
with our exact two-loop coefficients, we find that, with iVj = 4, the two-loop term in 
m/m(m) is overestimated by only 9% in ( P-1C| ), whilst in ( |3.12| ) one obtains an underesti- 
mate by 31%, and in ( |3.20|) an overestimate by 27%. We thus see some merit in estimating 
radiative corrections by naive non-abelianization, in cases where a large- Nf calculation 
is practicable, whilst an exact one is not. It is noteworthy that naive non-abelianization 



of the large- terms |26| in e + e annihilation underestimates the three-loop terms in R 



by less than 6%, with 4 or 5 active quark flavours. 



4 All— order results 



We now apply the methods of |3(| |3T|| to study the matching coefficients, at all orders 
of perturbation theory, in a highly fictitious limit: Nf — > — oo. 



As in other recent studies [15|, [u], [Ty|| , the intention is to investigate so-called renor- 
malon effects, associated with factorial growth of the coefficients of perturbation se- 
ries 



32|j . Unfortunately, there is no gauge-invariant prescription that enables one to 
obtain all-order results in the asymptotically-free theory of interest; instead one merely 
studies a single chain of fermion loops, such as that which gives rise to the Landau pole in 
QED. The hope is that one may learn something about the non-abelian theory, by imag- 
ining that asymptotic freedom still holds with an infinite number of massless flavours. In 
practice, the analysis amounts to no more than changing the sign of Nf in the large- Ay 
methods of @ [HJ, [Hfl. 



4.1 Master formula 

Following recent practice, we disguise the sleight-of-hand, by hiding the large- Nf contri- 
bution to the /5-function in a (fictitiously) positive value of (3q = |(11Ca — ATpNf). In 
other words, we imagine that (3 = — dloga s /dlog/i = f3oa s /27i + 0(o%) remains finite and 
positive, so that the QCD behaviour (3 ~ l/log(/x/Aj^g) still holds at large fx. 

In the matching coefficient we retain only the leading terms, of order l//3 , in the 
limit (3q — > oo, whilst retaining all powers of (3 ~ l/log(/i/Aj^g). As the heavy-quark 
loop is negligible, in comparison with loops from a large number of massless quarks, we 
immediately obtain r = T B and Zq — 1. A single L-loop diagram contributes to r . 
It is obtained by inserting a chain of L — 1 massless-quark loops in the gluon line of 
diagram |l|a. A corresponding insertion yields the L-loop diagram contributing to Zq. 
As (3q — > oo, the multiplications in ( |1.4| ) degenerate to mere addition and subtraction 
of terms of order l//3 - Hence multiplicative renormalization of the QCD and HQET 
currents amounts to no more than minimal subtraction of powers of 1/e. Moreover, there 
is no mass counterterm to consider at 0(1/ (3q), and coupling-constant renormalization 
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amounts only to f3 g 2 /(4Tr) 2 = fi 2e /3/(2 + f3/e), with /i 2 = // 2 e 7 /47r. Hence the entire 
perturbation series may be written, formally, as follows 

C r (ji) = 1 + £ F<k£ ^ £ " > (^ 2~f^ j ~ (minimal subtractions) + O (j^J , (4.1) 

where F(e, u) is regular at e = u = 0, and is easily calculated, in terms of a massless 
one-loop integral, in the gluon proper self-energy, and two massive on-shell one-loop 
integrals, determining the contributions of Zq and r . 



As in the case of deep-inelastic sum rules at large Nf pq| , the computational burden 
is slight, in comparison with the analysis of F32 hypergeometric functions in the two-loop 
integrals required for current correlators in QED |3(J and HQET |15|]. We find that 



^ C F f f i\^ T(l + u)T(l-2u) N(e,u) 

p Km) e -7*r(3-u-e) [DCe)] 1 -/*' 1 ] 

N(e, u) = (3 - 2e)(l - u){l + u - e) + 2 - u - e + 277(71 - 2 + e)u - 2(n - 2 + e) 2 ,(4.3) 
D( e ) = 6e 7£ r(l + e)B(2 - e, 2 - e) = 1 + \e + 0(e 2 ) . (4.4) 

The first term in ( |4.3|) derives from Zq; the remaining terms, from To, are quadratic 
in h = r](n — 2 + a), since they result from an essentially one-loop calculation. The 
function (|4.4| ) comes from the gluon self-energy. 

We now follow the methods of |29], |3(|, expanding F(e, Le) in powers of e and Le, and 
expanding [/3/(2e + (3)] L in powers of (3/2e, to obtain a quadruple sum in (|4.1|). As shown 
in |p9| , combinatoric identities relate \ je terms, and hence MS subtractions, to the Taylor 



coefficients of F(e, 0). In the case of the matching coefficient (|1.4|), we simply obtain 



7j - U = PF{-0/2, 0) + O ( . (4.5) 



As shown in ||30|1 , the finite terms receive contributions from the Taylor coefficients of 
F(e, 0) and also from the Taylor coefficients of -F(0, u). The former are scheme-dependent, 
and give a well-behaved series in (3. The latter are scheme-independent, and give a series 
that is not Borel-summable. A formal statement of the result may be written very simply: 

F(0,0)-FM) , J A f-2u\F(0,u)-F(0,0) , _ /l\ 



C r (/i) = 1 + / de + / duexp — — v ' ' + O 



(3 u \(3 { 



(4.6) 

The function F(0, u) specifies the Borel transform of the scheme-independent con- 
tributions flT5| , |TB|. Attempting to undo the Borel transform, via the notional Laplace 
transform of ( f4.6|) , one encounters singularities at u > 0, which are referred to as infrared 
renormalons [p2| . They reflect the fact that at higher and higher orders in perturbation 
theory one is probing regions of smaller and smaller gluon momenta. The first infrared 
renormalon occurs at u = ~, giving a singularity in the integral with a residue that is 
a multiple of exp(— l//3)/i/m ~ hrg^/m. Any attempt to sum the perturbation series 
involves an arbitrary choice of how to deal with this singularity. In general, a singularity 
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at u = uo > has a residue that is a multiple of (Aj^/m) 2u °. Thus the perturba- 
tion series, taken to all orders, has ambiguities that are formally commensurate with the 
higher-dimension operators in the 1/m expansion (|1 . 1| ) . 

The existence of such ambiguities is profoundly unsurprising. We see from (|3.13|) 
that the introduction of an infrared regulator, A, such as a gluon mass, would modify 
the result of on-shell integrals at order X/m. The infrared renormalon at u — \ serves 
to remind one of this simple fact. The situation is in close analogy with the operator- 
product expansion of QCD current correlators, at large Q 2 = —q 2 , where the first infrared 
renormalon, at u — 2, leads to an ambiguity of order (Aj^g/Q) 4 , commensurate with the 
gluon-condensate contribution [32]. Thus one should not flinch at renormalons; they serve 
as healthy reminders that one has chosen to integrate over all gluon momenta, including 
the infrared region where non-perturbative effects are dominant. In fact, infrared renor- 
malons have a positive virtue: the pattern of their residues must match the contributions 
of higher-dimension operators and thus provides consistency checks on the form of the 
1/m expansion of HQET, or the 1/Q 2 expansion of massless QCD. Since we have not 
taken the (inordinate) trouble of using an infrared regulator, the factorization of short- 
and long-distance physics into coefficients and operators is intrinsically ambiguous; each 
resummation prescription corresponds to a different set of values for the operator ma- 
trix elements. Nature, however, is not as slip-shod as we; she takes care to ensure that 
physical quantities are independent of the exigencies of our calculations. 

Hence we see, from ( |4.5| ) and ( |4.6| ), that the simple multinomial ( |4.3| ) furnishes de- 
tailed information about anomalous dimensions, through its e-dependence, and about 
renormalon singularities, through its ^-dependence. Moreover, it does so for all current 
matchings, with n specifying the QCD current and rj = ±1 the specific components that 
are matched to HQET currents. We now unpack some of this wealth of information. 



4.2 Anomalous dimensions 



Since we know that 7j vanishes at n — 1, we can obtain both jj and 7j from 



(n- l)(d- 1 -n) 



2tt 185(d/2, d/2)B(d/2 + 1, 3 - d/2) 



1j 



d=i+j3 



n=0 



(4.7) 



at order l//?o- For the QCD currents, we verify the n = result of [£9|, and the n = 3 
analysis of p6[| . The n = 2 and n = 4 results, as well as that for jj, are, we believe, 
new. It is rather intriguing that jj — -7^ vanishes at O(l//5o), thereby continuing, at all 
orders, a trend already apparent in the simple iV/- independent two-loop formula given 
for this combination in 0. A consequence of the all-order result is that the coefficient 
of the light-quark condensate, in the operator-product expansion of the HQET current 
correlator, is scale-independent, at leading order in 1/ ' (3q. This is confirmed, at two loops, 
by |[27|| . The all-order result may be somewhat accidental, since all large- Nf anomalous 
dimensions are fairly simply related. For example, we find, by similar methods, that the 
QCD field, q, and the HQET field, Q, have the following large-A^y off-shell anomalous 
dimensions: 

lg = a C F ^-i(3jj, ^ Q = aC F ^+ l 1 (4 + (3^ J , (4.8) 



7T 
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at order l//3o- The dependence on the gauge parameter, a, is limited to the one-loop 
level. In the Landau gauge, with a = 0, the light-quark result agrees with |33j and 



the heavy-quark result with | I5| . Impressively, all the 0(1//3q) terms in the anomalous 



dimension of the electron field were obtained in the QED analysis of [ 551 . 

We have a strong check of ( |4.8|) . The HQET current anomalous dimension is 

7J = 7r + |(7<? + 7Q), (4-9) 

where 7r = dlogZr/dlog/i, Zr being the renormalization constant of the HQET proper- 
vertex function. It can be obtained from the 1/e pole of the proper vertex at zero light- 
quark momentum and non-zero heavy-quark residual energy, because there are no in- 
frared divergences in this case. At order l//?o, all vertex diagrams with more than one 
loop contain a transverse gluon propagator, with quark-loop insertions. The argument in 
Section 2.1 shows that such diagrams vanish. Therefore, at O(l//3o), there are no multi- 
loop contributions to 7r- Moreover, the one-loop contribution vanishes in the Landau 
gauge. Taking account of this, the field anomalous dimensions (|4.8j) reproduce jj. 



One easily obtains the perturbative expansions of (|4.7|) and (|4.8j ), from the e-expansion 



v^oo ( 1 2k\ k 

2^k=0 \ 2 k 3 J t 



18B(d/2, d/2)B(d/2 + 1, 3 - d/2) exp {£^ 3 (3 + [-l] s - 2 s ) C(s)e s /s} 
= 1 - \e - if. 2 + (2C(3) - f ) e 3 + ( 3C (4) - f ) + O 



(4.10) 



which reproduces the large-iVy features of all existing two- and three-loop calculations, 
and provides checks for future perturbative calculations. It is notable that C(4), which is 
conspicuously absent from three- and four-loop QCD results, is bound to occur at higher 
orders. Whilst the perturbative expansion ( |4.10| ) is valid only for |e| < |, corresponding 
to \/3 1 < 1, one can see from the singularities of the Euler-Beta functions in ( |4.10|) that 
in fact no singularity is encountered for d > — 1. Hence the anomalous dimensions are 
well defined for (3 > —5, and summation of the large-iV/ perturbation series results in an 
extension of the domain of convergence. 



4.3 Renormalon ambiguities and their cancellation 



In stark contrast to the first integral in ( f4.6|) , which exists for all (3 > —5, the second 
contains infrared renormalons, for all > 0. If one arbitrarily chose to evaluate it by 
principal-value integration, one would differ, at order A^/m, with someone who chose 
to add some (presumably real) multiple of the residue at u = \. A solution would be to 
absorb this disagreement into different values adopted for the other terms that occur at 
0(l/m) in QLl) . 

The 1/m corrections to the ground-state meson matrix elements were calculated in [H| 
for the currents with n = 1. We performed the corresponding calculation for an arbitrary 
T, obtaining a general result in terms of h\ e=Q = i](n — 2): 

^4r^r = c ^ + 1 ( Gl + 2d ^ ~ - -)) + o (-1) , (4.ii) 

<0|J(/i)|M> m \ 6 / \m 2 J 

d r = l- 2r](n - 2) , d T = \ (cP r - 3 
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where M is the ground-state or excited-state meson that couples to the current, and G\ 
and G2 are bilocal matrix elements of the heavy-quark kinetic energy and chromomagnetic 
interaction, respectively. Multiplying T by 70 changes the sign of dr, without changing 
d r , since the latter is obtained from a^Y^a^^ = 2d T T^. We find that: d T = 
(— l) n+1 dr for all 8 currents; d-p = 3 for the 4 currents that couple to O 1 * 1 mesons; d-p = — 1 
for those that couple to 1^ mesons. The ground-state S-wave 0~ and 1~ mesons remain 
degenerate at this order in 1/m, as do the excited-state P-wave + and 1 + mesons. 

Viewed from the standpoint of HQET, the residual-energy term, A = m-Q — m^ + 
0(l/ra&), and the matrix elements, Gi } %, suffer from ultraviolet-renormalon ambigui- 
ties ||15|| . HQET knows nothing about the pole mass, m&, since it deals only with residual 
energies. However, the HQET self-energy, £(a>), has a linear ultraviolet divergence, con- 
ventionally suppressed by dimensional regularization. If one were to use an ultraviolet 
momentum-space cut-off, it would be necessary to introduce a residual mass, into which 
this linear divergence could be absorbed. Moreover, G li2 are bound to acquire ultraviolet- 



renormalon ambiguities, via mixing with A [17|. We now determine these by demanding 
that they cancel the infrared-renormalon ambiguities in (|4.11 ). 



Setting e = in the master result 



and 



we obtain 



F(0,u) = — f 




B(l + u,l-2u) 
2-u 



5 - u - 3u 2 + 277(71 - 2)u - 2(n - 2) 2 



(4.12) 

Comparing the residue at u = ~ with that for the pole mass |j!5 , 1^] , we obtain the 

ci 

of the pole-mass ambiguity, Am: 



ambiguity of the matching coefficients due to the infrared renormalon at u = |, in terms 



AC r (/i) 



+ 77(71 -2) -2(n-2) 2 



Am 



m 



(4.13) 



at any scale, \i. This ambiguity must be compensated in ( [4.1 1| ) by ultraviolet-renormalon 
ambiguities AGi and AG 2 Since ( 4.13|) is quadratic in 77(71 — 2), we have 3 equations, 
with only 2 unknowns. They are indeed consistent, and yield 



AG 1 



I Am , 

4 



AG, 



lAm. 



(4.14) 



This result was obtained in |I7[ , where the matching coefficients of the 77 = ±1 components 
of the n — 1 current were considered. However, the 2 ambiguities ( |4.14|) were obtained 
in JT7] from only 2 equations, with no consistency check. We have one check resulting 
from one extra equation. 

Finally, one might wonder what are the prospects of going to 0(l//3 2 ), where the 
renormalon structure is presumably far less trivial, with the possibility of cuts appearing 
in the Borel transform, rather than mere poles. The prospects appear to be rather slim, 
since one would have to insert chains of light-quark loops into the gluon lines of the 
diagrams of Fig. [1]. Even worse, one would need to insert a loop into the three-gluon 
vertices. There is, however, one term that is clearly tractable: that obtained by replacing 
the one-fermion-loop boson self-energy of diagram |l|b by self-energy terms with L loops, 
of which L — 2 are fermion loops. This was achieved in an 0(1/Nj) analysis of the muon 
anomaly [R(J, where the large-L behaviour was obtained, by evaluating an integral of 
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the Borel transform of the corresponding contributions to the Gell-Mann — Low function, 
also obtained in closed form in |3(J. Recently, corresponding contributions to the HQET 
self-energy have been analyzed |34| , with a result identical to that in |30 . 



5 Summary 



Our principal result is formula ( |2.19 ), which matches QCD and HQET currents in ( |1.3| ). 



It applies to all currents with a Dirac structure T = 7^ . . . 7^™' that anticommutes or 
commutes with 70. The one-loop term agrees with [|I|. The coefficients in (|2.20|) and ( 2.2 1|) 



give our new two-loop term, for any gauge theory, current, and quark-mass ratios. To 
use the general result one sets n to the number of 7-matrices in T and uses rj = +(— l) n , 
or 7] = — (— l) n , according as whether T anticommutes or commutes with 70. Multiplying 
T by a naively anticommuting 75 does not change the matching coefficient. Table 1 gives 
specific numerical results, ignoring small effects of light/heavy quark-mass ratios. 

All-order results, of similar generality, were obtained from the master formula ( |4.2|) . 
which gives the O(l//?o) terms in the matching coefficients ( [4.6[ ) and the anomalous di- 
mensions ( [4.7|) . Infrared-renormalon ambiguities in the matching coefficients are cancelled 
by the ultraviolet-renormalon ambiguities (|4.14j ), derived in |l7j without benefit of our 



consistency check. 

We briefly note the following salient points. 

1. Consistent results were obtained by including heavy-quark-loop effects in all 6 of 
the terms of (|1.4|); one must not omit heavy loops from HQET. 



2. The HQET on-shell renormalization constant ( |2.4|) removes the non-uniformity, at 



zero light-quark mass, of the corresponding QCD result |L2|, as shown in 



3. The current-independent procedure of ( |2.8| ) and ( |2.9| ) is an efficient alternative to 
evaluating traces with complicated antisymmetrized products of 7-matrices. 

4. Effects of non-trivial quark-mass ratios are obtained from the 3 integrals defined 
and evaluated in Appendix A, where series expansions are given. 

5. Apart from such two-scale integrals, our calculations were purely algebraic, as ex- 
emplified by the polynomial coefficients of Appendix B for the QCD vertex ( |2.13|) . 

6. All two-loop anomalous dimensions of QCD currents are given by ( p,17|) , whose 
specialization ( |2.18| ) to the tensor current appears to be new. 

7. The matching coefficients (|2.19|) confirm the universality of the finite renormaliza- 
tions ( |3.3| ) and ( |3.5| ) that restore chiral symmetry to the 't Hooft-Veltman 75 [24 . 



Two-loop corrections to ratios of meson decay constants are given: in ( |3.10|) , which 



confirms the relation between MS and pole masses |T]J; in fl3.12| ), which gives a 
— (3.1 ± 0.5)% two-loop correction to /b*//b, comparable to the —(4.6 ± 0.4)% 
one-loop correction; and in (|3.20|) , which gives the tensor coupling of B*. 



20 



9. 



All our two-loop corrections are distressingly large; fastest apparent convergence of 
/b* / /b would require evaluating a s at // = 370 MeV. 



10. Infrared safety of ratios of matching coefficients, at finite orders of perturbation 
theory, is exemplified by ( |3.13| ) , which also exposes the infrared renormalon at u — \ . 

1 1 . Decoupling of the t-quark loop from /b* / /b is explicitly demonstrated in (|3.14j) . 

12. The fe-quark loop approximately decouples from /b*//b, as shown by fl3.17| ). 

13. Naive non-abelianization of the massless-quark-loop contributions, by the process 
Ni — > Ni — y, approximates our exact two-loop terms, in ratios of decay constants, 
at the 30% level, or better. 

14. In the MS scheme, every large-iVy series of the type ( |4.1|) is formally resummed 
by Q4.6D , whose first integral gives the scale-dependence, whilst the second contains 
scheme-independent renormalon singularities. 

15. Anomalous dimensions and renormalon residues, for any current, are encoded by 
the e- and ^-dependencies of the master multinomial 



16. The anomalous current dimensions (|4.7|) involve the universal [T5|, |26], ^9| gOj, ^3 
^-expansion ( |4.10j ). 

17. The anomalous field dimensions (|4.8|) confirm our result 7j = \^q q + 0(1/0%). 

18. The 0(l/m) contributions for any current matching are given by (|4.11| ). 

19. The residue in ( |4.12|) at u — ^ produces the ambiguity (|4.13|) , which is cancelled 
in ( [4. 11| ) by ambiguities in the pole mass, m, and the matrix elements Gi.2- 

20. The inclusion of 0(1/(3%) terms in all-order calculations will be very demanding, 
though some progress in the abelian case has been made in [BUI E^, [34)1. 



In conclusion: we hope that our results will serve both to enable more accurate extrac- 
tion of meson decay constants from lattice simulations and sum rules, and also to provide 
detailed testing grounds for perturbative and non-perturbative approximations based on 
large-iV/ expansions. 
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A Dilogarithmic integrals A&(r) 



The finite-mass corrections (|2.6| ) and (|2.16|) involve 3 integrals of the polarization operator 
subtracted at m ; = 0. With z = —m 2 /k 2 , this takes the form (Tl[ 



U(z) = 2(1-2 z)Vl + 4z arccoth y/1 + 4z + log z + , (A.l) 



which is then integrated over simple functions ofy = 2/(l + -v/l — 4m 2 /k 2 ), to obtain 



i 



A 1 (r) = -J dy^—V-Il [r 2 ^r-\ = "(1 + r)L + (r) - (1 - r)L_(r) + log 2 r + C(2), 



1 f , 4-y / 2 l-y 

i 



-r(l - r 2 )L + (r) + r(l - r 2 )L_(r) + 2r 2 (logr + 1), (A.2) 



o 

A 3 (r) = dyyll [r 2 ^—^-\ = -r 3 (l +r)L + (r) +r 3 (l -r)L_(r) - r 2 (Yogr + ^ 



in terms of the dilogarithmic integrals [|11], |12 



o 



| loga:-logr _ J | log 2 r - log r log(l ± r) + ^C(2) - Li 2 (^r), r<l, 
' a;±r ~ \ logr log ^ + Li 2 ( T l/r), " r > 1 , 

(A.3) 

where Li p (x) = Y^=\X n /n p . Expanding ( |A.2| ), for r < 1, we obtain 

oo oo 

A fc (r) = 2 logr £ g k {n)r 2n + £ <^(n)r 2 " + 5 k , (A.4) 

n=l n=l 

*(») = 2^T-^. * = 3 « 2 ^ 

^ = 2^1-2^3' ^-3 C (2 M l-r 2 ), 

^ s(n) = 2^3 " 2^4 (n ^ 2) ' ^ = 1 ' ^(2) = -2, 
5 3 = - r 4 log 2 r + 3C(2)r 3 (l - |r) , 

where <7&( n ) = dgfc(n)/dn. Similarly, for r > 1, we find 

OO OO 1 

A fc (r) = -2 logr ]T &(n)-^ + £ 0k(«)-^ + > ( A - 5 ) 

n=0 r n=0 r 

^ l(n) = ^"^TT » 1 (°) = - 1 ' ^(0)=2, 5 1 =log 2 r + C(2), 

^ ) = 2^-2^TT' * 2 = ' 

^ ) = 2^T4-2^' h = °- 
Finally, external-flavour contributions, with = m, are obtained from 

A 1 (1) = 2C(2), A 2 (l) = 2, A 3 (1) = C(2)-|. (A.6) 
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B QCD proper vertex 



Here we present all non-zero coefficients in the formula ( |2.13| ) for the QCD proper vertex. 
The DTj^X source of these formulae has been generated by a REDUCE program, using 
the library package RLFI by R. Liska ]20 . 

[3d 3 - 30d 2 + lOld - HO) (d - 1) (d - 2) (d - 5) (d - 6) 
-2 (d 5 - 17d 4 + 109d 3 - 325d 2 + 428d - 172) (d - 5) (d - 6) 



«000 

a 001 

«002 

^003 
«010 

&011 
«012 
^013 
a 030 
«031 
«032 
&100 
a 101 
Ql02 
a 103 
«110 

am 

a 112 
«113 
(*120 
«121 
^122 
^200 

«201 
^202 
a 230 
&231 
a 232 



_ 

-2 (5d 4 



50d 3 



164d 2 - 156rf - 56) (d - 5) (d - 6) 



(d-4) 2 (d-5) (d-6) 

d 2 - 8d + 10) (2d - 7) (d - 2) (d - 4) (d - 5) (d - 6) 
42d - 32) (d - 4) 2 (d - 5) (d - 6) 



-2 (d 3 - 
-2 f 3d 2 



13d 2 

- lOd - 4) (d - 4) 2 (d - 5) (d - 6) 
8(d-4) 2 (d-5) (d-6) 

-2 (7d 3 - 76d 2 + 255d - 246) (d - 2) (d - 3) (d - 4) 
4 (3d 4 - 50d 3 + 287d 2 - 660d + 492) (d - 3) (d - 4) 
32 (d 2 - 6d + 6) (d - 3) (d - 4) 2 

- (d 4 - 9d 3 + 26d 2 - 24d + 2) (3d - 8) (d + 2) 

4 (2d 6 - 31d 5 + 183d 4 - 508d 3 + 642d 2 - 260d - 32) 

4 (2d 5 - 2d 4 - 139d 3 + 794d 2 - 1600d + 1088) 
-16 (2d - 7) (d- 2) (d-4) 

- (2d 4 - 24d 3 + 105d 2 - 199d + 134) (3d - 8) (d - 2) 

2 (2d 5 - 33d 4 + 210d 3 - 643d 2 + 916d - 448) (d - 4) 

2 (8d 4 - 71d 3 + 211d 2 - 178d - 88) (d - 4) 

-8 (2d - 7) (d - 3) (d - 4) 

4 (3d - 8) (d - 2) (d - 3) (d - 4) 

-16 (d- 2) (d-3) (d-4) 2 

-16 (d- 2) (d-3) (d-4) 

- (d 3 - 12d 2 + 45d - 46) (3d - 8) (d - 6) 
2 (d 3 - 18d 2 + lOld - 184) (d - 2) (d - 6) 
8 (d 3 - 12d 2 + 47d - 56) (d - 6) 

2 (d 2 - 5d + 2) (3d - 8) (d - 4) 
-4 (d 3 - lid 2 + 30d - 16) (d - 4) 
-16 (d- 2) (d-4) 2 
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